KOBE-TH-12-04 
HRI-P-12-09-001 
RECAPP-HRI-2012-010 



Quark mass hierarchy and mixing via geometry 
of extra dimension with point interactions 

Yukihiro Fujimotoj^ Tomoaki NagasawaJ^ 
Kenji Nishiwakij^ and Makoto Sakamoto^ 

^''^ Department of Physics, Kobe University, Kobe 657-8501, Japan 

^Tomakomai National College of Technology, 443 Nishikioka, Tomakomai 

059-1275, Japan 

^Regional Centre for Accelerator-based Particle Physics, 
Harish- Chandra Research Institute, Allahabad 211 019, India 



Abstract 

We propose a new model which can naturahy explain origins of fermion generations, quark 
mass hierarchy, and Cabibbo-Kobayashi-Maskawa matrix simultaneously from geometry of 
an extra dimension. We take the extra dimension to be an interval with point interactions, 
which are additional boundary points in the bulk space of the interval. Because of the 
Dirichlet boundary condition for fermion at the positions of point interactions, profiles of 
chiral fermion zero modes are split and localized, and then we can realize three generations 
from each five-dimensional Dirac fermion. Our model allows fermion flavor mixing but the 
form of non-diagonal elements of fermion mass matrices is found to be severely restricted 
due to geometry of the extra dimension. The Robin boundary condition for a scalar leads 
to an extra coordinate-dependent vacuum expectation value, which can naturally explain the 
fermion mass hierarchy. 
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1 Introduction 



Recently the ATLAS and CMS experimental groups of the CERN Large Hadron Collider (LHC) 
have announced the excess at 125 GeV, which is consistent with the Standard Model (SM) Higgs 
boson, with a local significance of 5a after combining 7 TeV and 8 TeV data [H [2] . This amazing 
thing means mysteries behind the last missing piece of the SM are ready to be unveiled. But the 
SM still has many points which are unclear in spite of lots of effort of physicists. 

One is called "quark mass hierarchy problem". In the SM, we are forced to comply with 
hierarchy of almost five orders of magnitude in Yukawa couplings of quarks for describing the suit- 
able quark masses. Closely related to this issue, the SM cannot answer the mechanism behind the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix, which describes the strengths of generation-changing 
interactions in the SM. In addition to the two issues, we cannot explain why we introduce three 
copies of quarks whose quantum numbers are the same except their masses and the degrees of mix- 
ings in the above interactions. Many attempts have been done for explaining the issues within four 
dimensional (4D) Quantum Field Theory (QFT) framework with, e.g., launching new continuous 
and/or discrete symmetry, introducing new matters and interactions, discussing renormalization 
group (RG) effects from a theory in a (very) high energy scale compared to the EW (Electro Weak) 
scale. 

When we focus on the case in five dimension (5D), where there is one additional spatial di- 
rection, we can find a new useful tool for tackling the above and the other problems, geometry. 
Two of the most renowned studies which show the power of geometry are [21 II], where the authors 
proposed innovative ways for solving hierarchy problem. Extra space can have a huge variety 
of structure, which are detected as differences in 4D effective theory point of view. In 5D QFT 
framework, we also find new mechanisms which we cannot find in 4D, e.g., generating spontaneous 
gauge symmetry breaking by global Wilson loop operator [3H7j, symmetry breaking by orbifold 
boundary conditions (BCs) (SHTU]. It, however, seems to be hard to answer the origins of the quark 
mass hierarchy, the quark mixing and the number of fermion generation only based on them. We 
comment that the existence of (compact) extra dimension(s) is suggested by the superstring theory. 
Lots of works have been done for settling the three problems in the quark (and lepton sectors) in- 
dependently and/or simultaneously in the many contexts of large extra dimension [HHIl], warped 
extra dimension P^HTTj . vortex profile [TSH2U] based on [21J, and others 

In this paper, we focus on one of the interesting mechanisms due to extra dimension, i.e., 
localization of fields. We can generate the hierarchy in the Yukawa coupling naturally when the SM 
fermions are localized at different points in one (more) extra dimension(s) ^25j, whose situation is 
realized by a 5D scalar field coupling to 5D fermions with a kink background [26l [27] . A variation 
of this possibility is to localize the Higgs scalar vacuum expectation value (VEV) in one extra 
dimension [28l 129] . Here we propose a simple way of realizing three chiral generations and their 
localization, where we introduce point interactions (or many branes) on an interval [30VI33] . This 
system is decomposed into multi intervals and due to the Dirichlet BC at the positions of point 
interactions, and fermion zero modes are split and degenerated. Each profile is localized around a 
corresponding point interaction as an effect of nonzero fermion bulk mass. 

When we construct a model with the above mechanism, it is reasonable to assume all the 5D 
fields live in the bulk with no tree-level localized term at the positions of point interactions. This 
setup is very similar to the minimal Universal Extra Dimension (mUED) model on S^/Z2 [M] . 
The mUED is one of the most investigated models in context of extra dimension and has many 
gripping points, e.g., existence of dark matter candidate is ensured by the accidental symmetry 
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under changing the two end points of /Z2 [35] |^ The latest concrete analysis for relic abundance 
of the candidate is found in Ref. ^Oj. The parameters of the niUED (and other six dimensional 
UED models) are restricted by the analyses based on recent LHC experimental results in Refs [H]- 
HS]. In the mlJED model, the Yukawa structure is the same with that of the SM, and therefore 
we still need the fine tuning in the coefficients. On the other hand, in our model, the 5D Yukawa 
couplings cannot possess generation index since one generation of the SM fields is introduced. In 
other words, another maneuver should be offered to overwhelm the difficulty. 

To generate the Yukawa coupling hierarchy via geometry, an extra coordinate-dependent and 
localized profile of the scalar VEV is preferred. An idea to realize this situation is to impose 
nontrivial BCs for the scalar field which is incompatible with its non-vanishing constant vacuum 
configuration. This mechanism has been applied to breaking translational invariance [IH], breaking 
supersymmetry [171 - I50] and extended to higher extra dimensions [5TI [52] . For the scalar singlet 
case, the profile is described with the Jacobi's elliptic function and we can find a parameter region 
where the elliptic function approaches the exponential function. The exponential form is ideal for 
generating large hierarchy within a natural choice of parameters and almost all the input param- 
eters take coefficients with (9(10) magnitude when we scale them based on each corresponding 
suitable mass value. 

This paper is organized as follows: In Section |2] we give a brief review on a way of constructing 
system with many point interactions briefiy and subsequently, debate about suitable choices of 
BCs for 5D fermion and vector fields. In Section |3] we search for the possibility for achieving quark 
mass hierarchy and mixing simultaneously in multi point interaction system with an exponential 
Higgs VEV profile. In Section |4] we construct a concrete model realizing the exponential VEV with 
high precision without violating gauge coupling universality and check the validity of the model 
through discussing naturalness of magnitudes of coefficients. Section |5] is devoted to conclusions 
and discussions. 

2 Basic properties of zero mode functions in system with 
point interactions 

2.1 Zero mode profile of 5D fermion on an interval with a point inter- 
action 

In this subsection, we consider the zero mode profile of 5D fermion on an interval with a point 
interaction, where is put in the middle of the whole system. A point interaction means that 
the interaction can occur only at a point as a 5-function potential. In Refs. [MHSS], possible 
point interactions in one- dimensional quantum mechanics are shown to be classified by boundary 
conditions which are characterized by U{2) parameters for a circle and f/(l) x f/(l) for an interval. 
According to this result, we will specify each point interaction by one of the possible boundary 
conditions in this paper. We use a coordinate y to indicate the position in the extra space and 
assign the locations of the three boundary points as (= Lq), Li, L2, respectively. The schematic 
diagram of our system in Fig. [T] helps our understanding. 

In this paper, we concentrate on a simple case, where there is no tree-level brane localized term. 

^ Recently a non-minimal version of UED model with brane-localized terms has been proposed in Ref. [36] and 
the collider physical studies on the model have been also done in Refs. [55Il39j . 
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Figure 1: A schematic diagram of the interval system with a point interaction aX y = Li. 



The 5D action of a fermion we consider is 



d X 



'-J0 



L2 



dy+ dy {idAiT^^ + Mp) ^ 



(1) 



where the latin indices run to 3, 5 (or y) and greek ones run to 3, respectively]^ and Mp are 
a 5D Dirac fermion and its bulk mass. We mention that 5D fermion bulk mass with the ordinary 
form is allowed in our system since we do not introduce orbifold Z2 parity. 

In what follows, we contemplate profiles of fermions. Due to variational principle, the following 
quantities must vanish at the corresponding boundaries as 



y=0 



y=Li-e 



0, 



0, 



y=Li+e 



(2) 
(3) 



where e is an infinitesimal positive constant. The form of \l/r^5\E' = can be decomposed into 
'^rS'^l = '^lS'^r = 0, where the 4D chirality is defined as 7^\Effl = ±\E'fl. Then, the Dirichlet 

L L 

BCs turn out to be consistent with Eqs. ([2]) and (|3]), i.e. 



R 



or = at ?/ = 0, Li ± e, L2. 



(4) 



We will, however, choose \E'ij = at all the boundaries (or \E'j;^ = at them) to realize the SM 
chiral fermions in the zero mode sector, as we will see later. 

We mention that once the BC of a right- (left-) handed part of a 5D fermion ^ is determined 

as 



^r{x,U) = Q (^i(x,L,) = 0), 



(5) 



where Li shows the position of a boundary point, the BC of the remaining left- (right-) handed 
part is simultaneously fixed through equation of motion (EOM) as 



{-dy + MF)^L = Q ((9, + M^)^R = 0) at y = 0, Li ± e, L2. 
We can expand the 5D fermion as 

■^{x,y) = ^!R{x,y) + ^!L{x,y) 



(6) 



(7) 



In this paper, we choose metric convention as tjmn = V 



MN 



matrices are T 



— T r = r?' 



-17 



diag(— 1, 1, 1, 1, 1). The representations of gamma 
^ ^ 'yO^i^^^a ^j-^j pg^y attention to the Chfford algebra is defined as 
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The series of {/^(n)} and {f^(n)} are eigenstates of the hermitian operators V^V and W\ respec- 

R L 

tively; the V (and V^) are defined as 



V:=dy + Mf 
:= -dy + Mf 



(8) 



(9) 



where M^i^^^ is a KK mass of the n-th right /left KK mode. This mass degeneracy is ensured by 
quantum mechanical supersymmetry (QMSUSY) [561459] as 



(10) 



For zero mode [M^,^^^ = 0), Eq. (10) takes the simple forms as 



Vf^,^{y) = Q 



(11) 



Taking account of the BCs, we find the zero mode solutions for \E'i = at y = 0, Li± e, L2 as (see 
Fig. 2) 



and 



4(0) (y) 



for = a.t y = 0, Li ± e, L2 as 

4(0) (y) 

^ 1 r 



and 



^ O T 








Mi e^'^y 









for < y < Li 

for Li < y < L2 

for <y < Li 

for Li < y < L2 ' 



for < y < Li 

for Li < y < L2 

for < ?/ < Li 

for Li < y < L2 



(12) 
(13) 

(14) 
(15) 



where Afi,Af2,N'l,J^2 normalization constants, whose concrete forms are shown soon later. 



Here we would like to comment on some important points of the situation. The Dirichlet 
boundary condition (|5| for fermion at y = Li turns out to effectively split the interval into two 
segments. Then, chiral zero modes are two-fold degenerate, and each profile of the zero modes is 
confined in one of the segments and localized exponentially at one of the edgesj^as shown in Eqs. 
(12) and (13) (or (14) and (15)). Thus, we can realize a model with two generations of 4D chiral 
fermions in the case of the interval with a point interaction. Furthermore, the localization of the 
zero mode profiles will be found to lead to the fermion mass hierarchy, as will see later. We can 
see pictures explaining this situation in Fig. [2j Each concrete form is written down as follows: 



^ It is noted that the value of Mp can take a negative value and in this case, the direction of wave function 
localization flips to the opposite side. 
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Figure 2: Profiles of tfie zero modes / (o) (y), (o) (y) (/. (o) (y), /. (o) (y)) are scliematically depicted 

^IR '^2R ^IL ^2L 

in (a) ((6)) for = = 0) with Mp > 0. Here there are the three boundary points at 
y = 0,Li,L2. 



• In the case of = at ?/ = 0, Li±e, L2, 

+ \/ e2M!^r_i ^''"^'""^^[^(^ - ^1)^(^2 - + (KK modes), (16) 

• In the case of = at ?/ = 0, Li±e, L2, 

+ /^^|S^^"''"^'^^^^[^(2/ - ^1)^(^2 - l/)]42(^)| + (KK modes), (17) 

where 6{y) is the step function and ipf''^ and ip^'' represent two (degenerated) fermion zero modes. 
ALj shows the length of the corresponding i-th segment which is defined as 

ALi = L,-Li_,. (18) 

Here every mode function is suitably normalized and we can find the factor for this purpose in 
front of the exponential functions. 
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We can also evaluate the right- or left-handed KK fermions profiles but the aim of this paper 
is to understand the mechanism explaining the fermion mass hierarchy. Therefore we pend this 
issue for a future work. 



2.2 Zero mode profile of 5D gauge boson on an interval with a point 
interaction 

Following the previous section, we move to the zero mode profile of 5D gauge boson on the interval 
with a point interaction. Here we concentrate on U{1) case since our intension is to investigate 
the structure of mass spectrum. The concrete form of the 5D action is 



dy + 



L2 



dy 



4' 



(19) 



where Fmn = d^iA^ — SatAm is the 5D field strength of the 5D f/(l) gauge boson Am and we 
assume that the background geometry is the same as in the fermion case. After taking variation, 
we focus on the forms of conditions at the boundary points as 



(dyA^ - d^Ay)SA>' 
{dyA^ - d^Ay)6A^ 



y=Li-e 



(dyA^ - d^Ay)SA^ 

J y=L2 

{dyA^ - d^Ay)6A^ 



0, 



0. 



y=Li+E 



(20) 
(21) 



Under the constraints, we could choose the following BCs, where the Neumann (Dirichlet) type 
BC is assigned for i^y), like the mlJED model as 



(dyAf,) = at y = 0, Li±6, L2, 
Ay = at y = 0, Li ± e, L2, 



(22) 
(23) 



in which A^'s zero mode with a constant profile is found. We mention that the compatibility 
between the BCs' (22) and (23) can also be shown from a viewpoint of QMSUSY [561 - I59] . But this 



configuration is problematic in a phenomenological point of view. In this setup, the zero mode 
profile of 4D gauge boson is also split ai y = Li and this fact means that a zero mode within an 
interval have limited interactions only with the particles inside the segment which the gauge boson 
belongs to (see Fig. 3). Since this possibility with doubly-degenerated zero modes {Ai^^ and ^42°^) 
is rejected in the SM, we have to change the BC aX y = Li. We further comment that the gauge 
universality in the SM is violated in this configuration and this fact gets to be another reason for 



discarding the system with the BCs given in (22) and (23). 



To avoid the problems, the profile would like to be continuous at y 
"continuous" conditions at this point for A^ and Ay as 



Li. Then we put the 



A, 



A, 



y=Li-e 



y=Li-e 



A, 



A 



y=L\+e 



y=Li+e 



and dyAfj^ 



and dyAy 



y=Li~e 
y=Li-e 



dyAfj, 



dyAy 



y=L\+e 



y=Li+e 



(24) 
(25) 



These conditions are consistent with the constraints in Eq. (21). In this case, we can observe only 



one zero mode of A^, whose situation is consistent with that of the SM (see Fig. [3]). It is noted that 
we need to put the continuity in the first derivative level because the Klein-Gordon (or Maxwell) 
equation is second order. We also consider this continuous condition for 5D scalar as in Eq. (pSj). 
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Figure 3: Profiles of zero mode 4D gauge boson witli BCs in Eq. (22) in (a) and witli BCs, wliicli 
are modified ai y = Li as in Eq. (24) in (b). 



2.3 Flavor mixing in system with point interactions 

We can also consider the "continuous" condition at point interactions for fermions as 



y=Li+e 



(26) 



where Li denotes a position of point interactions. An interesting application is given as follows: 
At this time, we add another 5D fermion ^' with the different bulk mass M'p and go to a two- 
point-interaction system, where the point interactions are located aX y = Li, L[ (see Fig. The 
5D action which now we think about is 



(27) 



where we do not divide the range of integral for y for clarity in description. The BCs for \E' and 
"if' are selected as in Fig. |4| where the red (blue) circular spots show the Dirichlet-type BC for 
left- (right-) handed part at the corresponding boundary points, respectively, and the small white 
dots mean the continuous condition for fermions. Under the BCs, zero modes of both \1/ and 
become two-fold degenerated and we distinguish the two states by adding new indices of "1" or 
"2" showing "generation" from the left to the right. When we choose the signs of the two bulk 
masses Mp, M'p as Mp > 0, M'p > 0, the way of localization of the zero modes is as in Fig. IS] The 



concrete forms of mode functions for \E' are the same in Eq. (17) and we know those for \E'' after 



adding prime symbol ' in some parameters as — )• i^'ilii^^L ~^ 2Li^f M'p,Li — > L'^ in 
Eq. (|T6|. 
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O 9 % % overview of system 

# # O # BC for ^ 

# 6 # # BC for 

— i i i i-^ y 

Li L\ L2 

Figure 4: This is an overview of our system with two point interactions, where the red (blue) cir- 
cular spots show the Dirichlet-type BC for left- (right-) handed part at the corresponding boundary 
points, respectively, and the small white dots mean the continuous condition for fermions. 




Figure 5: Zero mode profiles of 5D fermions \E' and ^' . 



The remarkable point in the system with the two point interactions is that the two zero modes 
with different 4D chirality and "generation" indices overlap in a middle region (Li < y < L[) 
of the system. This means that flavor mixing can be realized in our configuration. When we 
consider a UED-type scenario, right- and left-handed fermions are supplied by different 5D fields 
of SU{2)w doublet and singlet, respectively and separately. Thereby we can install this mechanism 
into UED model with no obstacle. In the SM, in fact, the flavor mixing structure is incorporated 
in the Yukawa sector with Yukawa couplings with the most general form and three copies of each 
SU{2)w doublet and singlet in the gauge eigenbases. These situations are affected exceedingly 
after adapting our mechanism for realizing flavor mixing. In addition, we should consider the 
profile of Higgs because it is a very important ingredient in Yukawa sector. These issues will be 
discussed concretely in the next section. 



3 Searching for the possibihty for achieving quark mass 
hierarchy and mixing in multi point interaction system 

Based on our discussion in Section [2| we make an attempt to create a model, where quark mass 
hierarchy and mixing are accomplished by use of the mechanism of flavor mixing due to multi 
point interactions. In this section, we consider the following form of action S with effective Yukawa 
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coupling in new geometry: 



Ls _ _ 

S = I d^x I dy{ Q {tdMT^ + Mq)Q+U {iduT^ + My) U 



Y'^''~^U{<P{y))U + Y^'^^D{(P{y))V + h.c. 



(2^ 



We have introduced one SU{2)w quark doublet Q, one up-type quark singlet U, one down-type 
quark singlet V with their 5D Dirac bulk masses Mg, M^, Mx>- We would like to emphasize that 
our model does not possess any generation index at this stage and that the fermion generation 
can appear dynamically. Here we have assumed that a "Higgs doublet" "H acquires VEV with 
y-position dependence such as 



H 







(29) 



and that the structure of the Yukawa sector is the same with that of the SM. Note that the 5D up 
(down) quark Yukawa coupling F*^"^ (y'^'^^) also do not contain any generation index of the quarks 
in our model. The SU{2)\y quark doublet can be decomposed as 



Q 



(30) 



Before we go into more details, we would like to point out some remarkable properties of our 
model. The form of the action in Eq. (28) seems to be very similar to the corresponding part of 
the mlJED model at the first glance. But there are two significant differences between the two 
theories, as explained below. 

One is about the structure of Yukawa coupling. In the UED, we should introduce three copies 
of fermion configurations to realize the three generation structure of the SM. Profiles of mode 
functions describing zero mode fermions take constant, therefore the fine-tuning in the Yukawa 
couplings is inevitable. On the other hand, we only introduce one copy of fermion configurations 
in our model. 

The other one can be seen in properties of the VEV of the Higgs boson. In mlJED models, BCs 
of the Higgs field are chosen as the Neumann-type, then the VEV gets to be a constant, whose 
situation is the same with the SM. Here we consider the Yukawa structure in our model briefly. 
As in Eq. (28), the 5D Yukawa couplings do not possess any generation indices. The SM Yukawa 



structure is expected to be produced through geometry of an effective form of Higgs VEV and the 
lopsided wave functions of zero mode fermions. In our model, the VEV profile of the Higgs scalar 
(0) is assumed to be ^/-position-dependent and to take the "warped" form of 



{(p{y)) = Aexp[a{y - L)], 



(31) 



with two massive parameters A and a, whose mass dimensions are 3/2 and 1, respectively. L (= L3) 
means the length of the total system. The reason for forming this conjecture is that this shape 
of the VEV is preferable to generate the large quark mass hierarchy in the SM within a natural- 
ordered setting of parameters. At this stage, we do not worry about the way of realizing this type 
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; L 



\H : 



for quark doublet (Q) 
for up quark singlet (U) 
for down quark singlet (V) 



Figure 6: A schematic diagram for explaining our attempt. Conventions of circular spots and 
dots are the same those in Fig. |4} 



of VEV profile and mainly devote our attention to phenomeno logical issues of it. In a later section, 
we discuss an example of materializing this (assumed) setup by use of the generalized Higgs BCs 
discussed in Ref. [5D] . 

We have to introduce new point interactions, where we take the Dirichlet BC, for all of Q,U,V 
to realize three generations in the SM. To this end, we assume that each 5D fermion feels nontriv- 
ially two point interactions with the Dirichlet BC. It turns out that zero modes of the fermions are 
three-fold degenerate, where each profile is confined in one of three segments. Here we consider 
the quark profiles explained in Fig. |6| where the meaning of the red and blue circular spots and 
the white dots are the same with those in Fig. |4} It should be emphasized that the positions of 
point interactions which each 5D fermion feels are not necessarily common. We assign the coor- 
dinates of the lower and upper ends of the total system with (= Lq) and L3, respectively, and 
the locations of point interactions between the two end points of the interval are represented by 
L^i\ L^2\ L^i\ L^2^ with the superscripts for identifying the type of the fields and the 

subscripts for showing the sequences when we count them from the left to the right (see Fig. |6]). 

Concretely speaking, we adopt the following BCs for Q,U,V, respectively as 

Qr = at y = 0, 4'^±£, L^'^ie, L3, (32) 
Ul = at y = 0, lJ"^ ± e, 4"^ ± e, L3, (33) 
Vl = at y = 0, L^f ± e, L^"^^ ± e, L3, (34) 

where three-fold generated left- (right-) handed zero modes emerge in Q {U, T>) as we have discussed 
before. Here we omit to describe the "continuous" condition for each field. It is noted that the 
profiles of U and D get to be the same in our configuration. The orders of the positions are settled 
on as 

(= Lo) < lS") < L'f^ < L^") < Li'^ < L3, 

(= Lo) < < < Lf < Li'^ < L3. (35) 

We comment that values of 4D effective Yukawa masses are evaluated as those of overlap 
integrals among the VEV, right- and left-handed modes. Then both magnitudes and signs of the 
bulk masses Mq, M^, M-p govern an important part of the results. 
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The forms of effective 4D Yukawa masses among the three generations are obtained after 
integration over y as 



Li 



dy 



Y^''^U{<P{y))U + Y'^'^^D{(t){y))V + h.c. 



diiix),d2i{x),d'-^l{x) 



U 



u. 



u 



(0)- 

IR 
(0) 
2R 
(0) 
3R 



(0) 
3iJ 



+ h.c., 
(36) 



where the mass matrices A^^") and A^^"') have the following structure 





" (n) 


{«) 







- (d) 

my 


(d) 
my 










{«) 
"^22 


(u) 







id) 

my 


(d) 
my 










(u) 

my 










(d) 
my 



(37) 



The three-fold degenerated zero modes are distinguished by the generation indices "1,2, 3" in both 
and down-type ones (d). Each matrix component in Eq. (37) is calculated 



up-type quarks (u) 

by overlap integrals among the Higgs VEV and zero mode functions of right- and left-handed 
fermions, which are functions of the fermion bulk masses Mq, Mu, Mx> and the locations of the point 

Llf'\ In contrast to the SM, some elements of the Yukawa mass 



interactions L^i \ Lg'"*, 



r H T (") T (d) 

-^1 5-^2 ^-^1 ! 



matrices are zero. The diagonal parts of ^22'', m^ ,m^ii ,m^^ ,m^^rj are always nonzero unless 



which, of course, is unsuitable for our 



we take an extremal parameter choice, e.g., Li 
purpose in this paper. We notice that the characteristic form of the mass matrices in Eq. (37) is 
given by the geometry of our setting. Some of the components in the mass matrices vanish due to 
no overlapping of mode functions. 

Which non-diagonal component is nonzero depends on the positions of the point interactions. 
Following the rule in Eq. (35), the only (1,2) and (2,3) elements are nonzero. When we reverse 
the order of position between h'f'' and L^"^ as 



(38) 



the value of the 



'1,2) component vanishes 



m 



(u) 
12 



0) but that of the (2, 1) component becomes 



nonzero (^2"^' 7^ 0). This issue is easily understandable when we focus on the fact that the 
row (column) index of the mass matrix corresponds to the generation of the left- (right-) handed 
fermion, respectively. We would like to mention that in our model, flavor mixing is inevitable to 



realize the SM quark masses. Because in the aligned configuration of L^i"* = L^^' = Lf , L 



-{u) 



<d) riq) 



L^2 \ it is very hard to realize the quark mass patterns of the first generation 



m 



up 



^ '^down) 



and those of the third generation (mtop > ?Tibottom) simultaneously. 

The point is whether we can reproduce the structure of the CKM matrix or not based on the 
limited forms of the Yukawa mass matrices, where we never find the overlapping between first 
and third generations. After some calculations, we can find a set of input parameters where we 
reproduce the quark mass hierarchy and the CKM matrix simultaneously. In the next subsection, 
we discuss, in detail, whether we can simultaneously reproduce the structure of the quark mass 



hierarchy and the CKM matrix based on the limited forms of the Yukawa mass matrices (37) 
without the overlapping between the first and third generations. 
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3.1 A solution in the setup with a warped VEV 

Let us first give the concrete forms of the Yukawa mass matrices in our system which have been 



discussed in the previous subsection. The fields Q,U,V with the BCs in Eqs (32), (33) and (34) 
are KK-decomposed as follows: 



Eti"!? (3;)/ (o)(?/) 



^(0). w'^ / X I +(KK modes), 



U{x,y) 
D{x,y) 

U{x,y) = + (KK modes) 

i=l 

3 

V{x,y) = Y^df^{x)f^o,{y) + (KK modes). 



(39) 
(40) 
(41) 



i=l 



where we only focus on the zero mode parts. Here the zero mode functions are obtained in the 
following forms 



f ,o)(y) = ui'^^e'^^^y-^^^^ e{y - L^^eiL^ - y) 

liL L 

/„,o,(y) = Ar(")e-*^"(^-^»'"'^) \9{y - Lt\)e{L^;^ - y) 

iR 

fMy) = Ml^^''-^-^^ \e{y - Lf\)e{Lf^ - y) 



(42) 
(43) 
(44) 



where we use the conventions, for clarity, of 



AL? 
(= Lo) = 
L3-- 

2Mq 



- (0 r (0 



Lili (for 2 = 1,2,3; / = q,u,d), 



-ii) 



-(<?) 



r(")_ rW 

rW_ rid) 
-^3 — ^3 5 



2M; 



u 



l-e 



2MuAL, 



(u) ' 



2M 



1 _ ^~2MTyAL\ 



(d) • 



(45) 



are wavefunction normalization factors for / (o) , / (o) , / ,(o) , respectively. The length 
of the total system L takes the universal value of 



T T _ t(<i) t{i) _ t{u) r(") _ rid) j(d) 

Ij .— 1^3 — ^3 — i^o — -^3 ~ ^0 — ^3 ~ ^0 ■ 

Here we choose the signs of the fermion bulk masses Mg, Mu-, M-p as 

Mq >0,Mu< 0, Mv > 0, 



(46) 



(47) 



where we take a negative value in Mu to generate a large overlapping for the top quark mass 



mainly) in . An outlines of the wavefunction profiles are summarized in Fig. 



12 



[up type] 



f (0) 



[down type] 




-(d) r (d) 

Figure 7: An outlines of the wavefunctions profiles. 



Each component of the A^*^"^ and Ai^"^^ is acquired by calculating the corresponding overlap 
integral as follows: 



mS^) = F(«) / rfyf (o)(l/)/„(o)(y)(0(y)) 



qL 



,{MQ + \Mu\+a)L'^ 



(u) 



MQ + \Mu\+a 



(4J 



(9) 92L "2iJ 



MQ + \Mu\+a 



(49) 



'(g) 



' - (9) VSL "Sfl 



g(MQ + |Mi,|+a)Lf _ ^{MQ + \Mu\+a)Li'''> 



(50) 
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(9) ('A,f„_Llj\/f,,l^„,Ar('") 



(9) 



I ^ (it) "Sfl 



2 3 ^ MQ + \Mu\ + a ^ > V ; 



-(d) 



rrig)=y('^) ldyf Mf^,.,{ymy)) 



2 2 ^ Mq- Mr, + a ' ' V ; 



(9) 



mg)=y('^) I dyf ^o^{y)f^,,{ymy)) 

' - (9) VSL "3ii 



r ^ 
^2 



^(J=FW \ dyf Mf^y^km) 

- (d) "2i? 



' - (d) «2L ''Sii 



^2 



where we have used the forms of wavefunctions in Eqs. (42), (43), (44), the conventions in Eq. (45) 



and the assumed profiles of the Higgs VEV in Eq. (31). In this paper, we choose the parameters 
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up quark 


mass 


down quark 


mass 


up (u) 
charm (c) 
top (t) 


1.8-3.0 MeV 
1.250-1.300 GeV 
172.1-174.9 GeV 


down (d) 
strange (s) 
bottom (6) 


4.5-5.5 MeV 
90-100 MeV 
4.15-4.21 GeV 



Table 1: Experimental values of quark masses from Ref. ^T\ . 



as 



Mo 
a 



0.338 • L, 
: 0.0115 -L, 
0.223 ■ L, 



t(i) 

r{«) 
^2 



8.67 -L 



-1 



0.689 ■ L, 
0.540 ■ L, 
0.676 ■ L, 
Mu = -7.98 • L~ 

y{«)/y(<i) = 12.0, 



L 



id) 



rig) 

^3 — 


1-L, 


Tin) _ 

^3 — 


1-L, 


rid) _ 
^3 - 


1-L, 


Mv = 


6.16 -L-i, 




= 275 GeV, 



(5^ 



which reproduce the quark mass hierarchy and the CKM matrix with a good precision. The quark 



mass eigenvalues and the CKM matrix are evaluated from Eqs.(48)-(57). The diagonalized Yukawa 
mass matrices take the forms of 



I diagonal = diag (2.13 MeV, 1.18 GeV, 174 GeV), 
A^^'^^ldiagonai = diag (3.85 MeV, llOMeV, 4.19GeV), 



and the CKM matrix is given as 



IK 



CKM I 



0.976 0.213 0.00448 
0.213 0.976 0.0475 
0.0145 0.0454 0.999 



(59) 
(60) 



(61) 



Subsequently, we present the latest experimental values. The quark masses are summarized in 
Table [l] and the CKM matrix is 



CKM exp. 



0.974 0.225 0.00415 
0.230 1.006 0.0409 
0.0084 0.0429 0.89 



(62) 



from Ref. pT]. From the results in Eqs. (59) and (60), the mass eigenvalues are well reproduced 
within about twenty percent differences and the absolute values of the CKM matrix elements are 
also described within about twenty percent differences in almost all the elements. But we should 
comment on the (3, 1) component of the CKM matrix which we have obtained, where the deviation 
ratio measures up to ~ 70%. 

At this stage, we ponder over the number of the input parameters and constraints originated 



from the equations. After keeping the relations about the positions in Eqs. (45) and (46) in mind, 
we count the independent d.o.f. (degrees of freedom) as 7 for the positions of boundary points, 
3 for fermion bulk masses, 2 for 5D Yukawa couplings, 2 for effective Higgs VEV parameters, 
respectively, and the total d.o.f. is 14. Meanwhile, the total number of the constraints is 9, 
where they consist of 6 quark mass eigenvalues and 3 CKM mixing angles. Here the two massive 
parameters of 



yW (or Y^^\ A), L 



(63) 
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are still not specified. It is noted that once one of Y^'^\Y^'^\A is fixed, the other is determined 
simultaneously. The reason for this indetermination is that the electroweak scale is not indicated 
within our analysis. 

The results in Eq. (58) are realized broadly within (9(10) range (when we choose the basis of 
scaling as L~^) and it is shown that we can explain both of the quark mass hierarchy and the 
structure of the CKM matrix in a natural statement. 



4 An example for realizing warped Higgs VEV by gener- 
alized Higgs boundary condition 

In the previous section, we have verified the possibility of accomplishing the quark mass hierarchy 
and the structure of the CKM matrix at the same instant without fine-tuning of parameters in the 
setup, whose geometry contains many point interactions. Here we have assumed the forms of 5D 
Yukawa couplings and the warped shape of the Higgs VEV profile {(f){y)). In this section, we give 
an example of generating this situation without conflicting with the physics at the SM. 

4.1 Generalized boundary condition for a scalar 

Following [20], we can consider the physics which is described by the 5D actions on a single interval 
[0, L] for a 5D scalar $ of 

S^ = j di/|<l>t9'^9/,$ + $t92$-y(|$|2)|, (64) 

with the bulk potential V^(|$p) = M^|$p + ||$|^ and where we do not consider tree level brane- 
localized term. The parameters M^, A need to be real due to hermiticity. We can find a difference 
from the ordinary UED-type form in the structure of the derivatives. Based on variational principle, 
the following form should vanish at the boundaries y = 0,L: 

^^dyS^ - {dy^)^^ = at y = 0,L. (65) 

It turns out that a larger class of BCs is allowed with two new real massive parameters L±, whose 
mass dimensions are —1, as 

$ + L+dy^ = at y = 0, 

$ - L_dy<^ = at y = L, (66) 

where L± can take values in the range of — oo < L± < oo. This type of BC is called Robin boundary 
condition. It is obvious that the above conditions include the ordinary Neumann (Dirichlet) BC 
in the case of L± = ±oo {L± = 0). 

4.2 Position-dependent scalar VEV in generalized boundary condition 



In this subsection, we discuss the profile of the scalar singlet $ and its phase structure. In the 
mlJED model, where the Higgs doublet takes the Neumann BC at the both boundaries, the 4D 
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effective Higgs potential is minimized with ease, where the VEV is a constant (y-independent) and 
the situation is the same with the SM. 

On the other hand in our model, the phase structure is nontrivial because the massive param- 



eters L± emerge in our generalized BC in Eq. (66). To investigate properties of phase structure of 



the scalar $, we have to solve the minimization problem of the functional indicating 4D effective 
Higgs potential 



^[$]:=^ dyj -<l>t52$ + r(|$n| = rfyj -$t52^ + M2|$|2 + (67) 

It is important to incorporate the kinetic term along the extra spacial direction into the potential 
because the minimum configuration of the singlet $ probably possesses y-dependence, as we will 
see below. 

Hereafter we search for the form of the VEV ($(y)) by solving the EOM 

{-d't + M2)<I> + A<l>"^$2 = 0, (68) 



which can be obtained after taking variation in Eq. (67). The solutions of Eq. (68) are found to be 
expressed in terms of the Jacobi's elliptic functions generally. The types of solutions with > 
are classified based on a parameter Q, which is an integration constant after integrating along y, 



with mass dimensions 5. It turns out that with the choice of Q < a solution of Eq. (68), which 
can realize a desired "warped" VEV in an asymptotic form, is give by 

my)) = ^ ■ . r- ^ ^' (69) 

tf(l/-l/o),/^^, 



cn 



where the several parameters are defined as 




^ 4A|Q|A 



(72) 



I/O means a translation d.o.f. along ?/, which appears after integration for solving the equation 



in (68). The index k is an important parameter of elliptic function to determine the profile. We 



can rewrite k by use of the input parameters in /i, v as 




where we conclude that the possible region of the value of k is 



2 < ^ < 1- (74) 
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Figure 8: These three plots represent the VEV profile of ($(?/)) = l/cn(?/, fc) in (a), its first 
derivative dy ($(?/)) in (b) and the form of ($(?/)) /dy ($(?/)) in (c). The red, blue, magenta, green 
curves are in = 0, 0.71, 0.9, 1, respectively. The red dots are the points oi y = 7r/2. 



Here we mention that the condition A > is required to ensuring the stability of the vacuum. We 
again write down the form of $ with less abbreviation: 



(*(?/)) 



M 

71 



jv^TTx- i} 



1/2 



1 



X 



cn (^M{l+X}^/^(y-yo),y[(] 



1 



(75) 



Next we decide to choose our strategy to treat scalar BC. We rewrite Eq. (66) as follows: 
= - ($(0)) /dy ($(0)) , L_ = + (<|.(L)) /dy (<|.(L)) . 



(76) 



The L^,L^ are determined by the values of ($(|/)) ,dy ($(?/)) at the endpoints y = 0,L and we 
discuss the shapes of ($(?/)) , dy ($(?/)) at first. In Fig. |8| we show three plots of ($(?/)) , dy ($(|/)) 
and {^{y))/dy ($(?/)), with suitable normalizations. The red, blue, magenta, green curves are in 
k = 0, 0.71, 0.9, 1, respectively and the red dots are the points of y = 7r/2. We give some comments 
in order. 



• The function of l/cn{y) can be represented as trigonometric or hyperbolic function in ex- 
tremal cases of k, concretely, which is l/cn(?/)|fc=o = 1/ cos{y) and l/cn{y)\k=i = cosh(y). 
Following the change of the value of k from to 1, the profile of l/cn(y) smoothly shift from 
1/ cos(?/) to cosh(y). 

• ($(y)) |fc=o and dy ($(?/)) |fc=o are divergent at y = 7r/2. Increasing the value of k from to 
1, this divergent point moves from 71/2 to infinity. 

• The profile of [($(?/)) /dy {^{y))]\k=o is divergent at y = tt and takes zero value a.t y = 7i/2. 
As increasing the value of k from to 1, these points move from tt (7r/2) to infinity. This 
profile is also divergent at y = independently of the value of k. 

• In the region between y = and y = y^, where y^^ is a first divergent point with positive value 
and, of course, ?/d = vr/2 in case of k = 0. The profile of ($(?/)) /dy ($(y)) is monotonically 
decreasing independently of the value of k. 

In this paper, we focus on the segment of {0,yd)- At this segment, as we review above, the profile 
of ($(|/)) /dy ($(?/)) is monotonically decreasing independently of the value of k and takes positive 
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IMI (L++L_] 



l+ArL+L_ 



unbroken phase 
<^>>= 



broken phase 







Figure 9: Phase diagram of the scalar $ in the case of L_|_L_ < 0, L_|_ + L_ < in Ref. |60j 



value. Therefore, the values of L+,L_ obey the following condition: 

(77) 




where the condition, L+ + L_ < 0, is automatically derived from the conditions of Eq. ([77j). The 
phase structure of the scalar singlet $ with the generalized BC is explored in Ref. [60] and we 
quote a phase diagram in our case of < 0, L+ + L_ < as Fig. [9j{} Lmax in Fig. [9] is defined 

as 

Lmax := max {L+,L_}= L_. (78) 

The important point is the existence of a critical length Lc in the region of > l/L"^^^, whose 
definition is 



Ml V l + M^L:L_ J LI 



max 



and when we introduce a U{1) gauge boson, the gauge symmetry is broken (unbroken) for L < Lc 
[L > Lc). On the other hand, the gauge symmetry is always broken in the region of < jj—- 
Now we discuss properties of the solution. On the condition of 

k^l and -y|^yo»l, (80) 
the form of the solution gets to be exponential as follows: 

my))^i^cosh(^fi{y-yoyj ^ ^g-v^''^^" ■ e^'^^ (81) 

which is just the form of the warped VEV. The detailed analysis including numerical calculation is 
executed in the subsection 4.4 For more concrete understanding, we go for another analysis, which 
is based on the discussion in Ref. [5D]. We introduce the eigenfunctions f\n){y) of the eigenvalue 
equation 

-dlhn){y) = E^n)fin){y), n = 0, 1, 2, ■ • • , (82) 
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with the EC's 

/(„)(0) + L+9,/(„)(0) = /(„)(L) - L^dyfUL) = 0. (83) 

In terms of the orthonormal eigenfunctions f{n){y), whose orthonormahty is ensured by the her- 
miticity of the operator {—dy), the field $ can be expanded as 

oo 
n=0 



with the corresponding coefficients Inserting it into £[^] in Eq. (67) leads to 



oo 



^[^] — I'^wl ~'~ (l^^rtic terms in 0(„)) , (85) 

n=0 

where 

m^„) :=M2 + E(,), n = 0,l,2,---. (86) 

Note that the quartic terms are non-negative for any configurations of (/)(„) because they come 
from the term of Jq dy^ |$|^ (> 0). It follows that the configuration of the VEV is given by 
($) = (or (0(n)) = for any n) if m^^^ > for any n. For realizing symmetry breaking, the 

condition m^-j < is mandatory and the |0(o)| term's contribution is probably dominant around 
the minimar] Consequently, we could approximate the form of ($(?/)) as 

(<l>(l/))~0(o)/(o)(l/). (87) 
Here we take the two parameters of the BC as 

7^ = T (M + e) , (88) 

where e is an microscopic (but not infinitesimal) and positive value. It is not difficult to obtain 
the form of /(o) (y) and the corresponding eigenvalue -E'(o) with the assumptions 

f(o)iy) ~ e(*^+^)^ E(o) = -(M + e)2, (89) 
where we derive the form of the warped VEV again. By calculating the value of m(o) as 

mjo) = -2eM - e^, (90) 

we can infer that the symmetry breaking is realized for a certainty with our premise of M > 00 

^ Continuing discussions in this way, we can classify the phase structure of the scalar without knowing the details 
of the VEV ($). 

^ In other words, the condition < l/L"^^^ is fulfilled. 



20 



4.3 A model for realizing our mechanism without violating gauge uni- 
versality in multi point interaction system 

Now we know that the form of the warped VEV can be achieved by the scalar with the generahzed 



BC in Eq. (66). Naively thinking, in the UED-type model with one Higgs SU{2)w doublet H with 
the generalized BCs at the two end point of the total system {y = 0,y = L), and the "continuous" 
conditions at the others, we can explain the quark mass hierarchy and the structure of the CKM 
matrix simultaneously via geometry. But two nontrivial issues exist in this setup. 

In the SM, the profile of the VEV {H) can be rotated by use of the SU{2)w global symmetry 

as 

^ r, /.%t1 . (91) 



V/V2L 

with V = 246 GeV because the profile (H) is constant. The first issue is that whether we can rotate 
the profile and obtain the above form or not. In our case the VEV becomes y-dependent generally 
and the 4D effective Higgs potential in Eq. (|67|) is considered to hold a very complicated structure. 



Thus it is nontrivial that whether the VEV {H) can take the form in Eq. (91). 

The second issue is very critical and important. When the VEV is y-position dependent, zero 
mode profile of a 4D gauge boson is determined by a Lame-type equation and is not constant any 
longer. Then the overlap integrals for quark- ant iquark-gauge boson interactions in the SM become 
generation-dependent and as a consequence, the gauge universality in the SM is jeopardized. Of 
course this result is unacceptable and thereby we need to alter our strategy. 

A remedy for this plague orders the coexistence of a Higgs doublet H and a singlet scalar $, 
where the BC of the former is chosen as the ordinary Neumann-type at the end points as 

dyH\y = = dyH\y = L = 0, (92) 



and of the latter is selected as the generalized BC in Eq. (66). The 5D actions for H and $ are 
given as 



Sh = I d^'x £ dyS^H^DMD^' + M")H - ^ {H^hY^ 
S^ = j d^x j\y^^^ {OmO^ _ m^) $ - ^ ($t$)' 



(93) 

(94) 

Here we introduce the 5D gauge bosons of Gm,Wm, Bm, which are 5D SU{3)c, SU{2)w, U{1)y 
ones, respectively. Dm means the covariant derivative for the corresponding gauge bosons. The 



5D fermions in Eq. (28) are also gauged for reproducing the SM interactions. Dm and the action 
for the 5D gauge bosons take the same form with those in the mUED and we do not write down 
them since the detailed information is not important in the following discussion. The BCs for 
Gm, Wm, Bm are chosen as 

dyGf^\y=o = dyG^\y=L = 0, (95) 

6*^1^=0 = Gy\y = L = 0, (96) 

where we only illustrate the gluon case. We write down a schematic diagram for explaining the 



BCs for bosons in Fig. 10 The green, orange and purple circular spots represent the ordinary 



Neumann, Dirichlet and the generalized BCs in Eq. (66), respectively. Except the both end points 
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Figure 10: A schematic diagram for explaining the BCs for bosons. The green, orange and purple 



circular spots represent the ordinary Neumann, Dirichlet and the generalized BCs in Eq. (66) 



respectively. Except the both end points of the total system, we put the continuous condition 



discussed in Eqs. (24) and (25). 



of the total system, we put the continuous condition discussed in Eqs. (24) and (25). Each 4D 



vector part has a zero mode, whose mode function is a constant. By putting the condition of 
M'^ > 0, the SU{2)w x U{1)y gauge symmetry is spontaneously broken through the usual Higgs 
mechanism and the SM situation is duplicated at the sectors of Gm, Wm, Bm, H- Here we should 
assign the U{1)y charge of the Higgs singlet $ as zero to make the situation that $ do not couple 
to any gauge bosons. Accordingly, the problem in gauge universality never occurs in the refined 
setup. 

What we should concern next is the structure of 5D Yukawa interactions. When we adopt the 
forms in the SM (or the mlJED), the large mass hierarchy cannot be created since the profile of 
{H) is constant. To simplify the situation, we introduce the following discrete symmetry 



(97) 



to prohibit the terms of Q{ia2H*)U, QHV, ^QQ, ^UU, ^VV with the Pauh matrix cts. The de- 
sirable 4D Yukawa structure is generated by introducing the following terms 



y 



d X / dy< $ 



y^''^Q{ia2H*)U - y^^^QHV + h.c 



(9^ 



where those operators are higher-dimensional compared to the previous five operators and allowed 
under the discrete symmetry in Eq. (97). We mention the coefficients 3^''"'', J^^'^-' have mass dimen- 
sion —2. After the ElectroWeak Symmetry Breaking (EWSB) occurs with nonvanishing {H) and 
($) , 4D effective Lagrangian which we assume in Eq. ( 28 ) is realized without any serious problem 
conflicting with the nature of the SM. 

Here we discuss some related issues. Whether 5D gauge invariance is conserved or not is one 
of important criterion for judging validity of BCs. When we break gauge symmetry by BCs, the 
issue of possible unitary violation due to longitudinal components of massive gauge bosons should 
be considered [621461] . We, however, mention that in the cases of the Dirichlet BC for the fermions 
at mid and end points and the generalized Higgs BC at end points, 5D gauge invariance is intact. 

In the above analysis, we neglect the contribution to the total scalar potential of the doublet- 
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singlet mixing term with a coefficient C 



^ ■ ■ 

'-'mixing 



(99) 



where the discrete symmetry in Eq. (97) cannot ban this form. After considering this part, the 



profiles of $ and H are deformed and the problem in gauge universality is revived. Therefore we 
should choose a sufficiently small coefficient C to avoid this obstacle. The detailed discussion is 
done in Appendix |Xj 

4.4 Detailed numerical calculations for justifying the model with Higgs 
doublet and singlet 

Based on the previous discussions, we reexamine the issue for validness of our model including 
the Higgs doublet and scalar singlet with numerical calculations. At ffist, we reconsider the 



approximation in Eq. (81). As we have discussed before, the cn function is almost equivalent to 



the cosh function in the limit of ~ 1, and then we obtain the form 



($(?/)) ~ u cosh 



- yo) 



ev 2 



t/'(j/-J/o) _|_ ^-y/^Kv-yo) 



(with k ^ 1). (100) 



If the condition y 2f^iy ~ Z/o) ^ 1 is fulffiled, the second term on the right-hand-side of Eq. (|100|) 
can be neglected within about 10% error (e~^ ~ 0.135). Therefore we get the outcome of 



my)) 



V 

~ — e 
2 



|m{j/-i/o) 



with /c ~ 1 and \\ -fiiy - yo) > I 



(101) 



Here we consider the situation of A; ~ 1 more concretely. As we show in Eq. (73), the parameter 



k is composed by some input parameters for the solution and /c ~ 1 is equal to the condition 



X 



4A|g| 

M4 



0. 



(102) 



It is obvious that for matching up to this condition, the smaller (greater) value of A and/or \Q\ 
(M) is preferred. But the extremal choices of A = 0, |Q| = 0, M = cxd turn into disorder 



and unnaturalness. We rewrite the approximated form in Eq. (101) with input parameters by 
considering the shapes of /i, u in Eqs. (70) and (71), which are roughly under the situation ~ 1 



'2 V^lQl 



(103) 



with the zeroth (ffist) order approximation in X for /i (z/). When we evaluate up to X's zeroth 
order, the value of u goes to zero and this is meaningless. By use of the results, we can rewrite the 



equation in (101) as 



my)) 



\Q\_]_ -Myo . My 

2 M 



(with /c ~ 1 and M{y — yo) > 1) . 



(104) 
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Figure 11: (a) the value of k in the function of X in the range of < X < 5000. (b) the value of 
the complete elliptic integral of the first kind K[k] in the range of 0.9 < A; < 1. 



The correspondence of Eq. (104) to the assumed warped VEV in Eq. (31) is as follows: 



A 



Ml 

2 



,M{L-yo) 



a 



M. 



(105) 



Then we notice that the shape of the approximated profile in Eq. (104) is mainly determined by 
M, I/O, and that A does not appear up to the approximation. 

In what follows, we discuss the validness of the above approximation and the deviation from 



it when we consider the exact form in Eq. (75). At first, we define the dimensionless parameters 
with tilde ~ in the basis of the massive parameter of the total length of the system L (= L3), e.g., 



M = ML- 



(9) 



(106) 



where part of the dimensionless parameters is already calculated in Eq. (58). The significant point 



is the bulk mass of the Higgs singlet is already almost fixed because of Eqs. (58), (104) and (105) 
as 



M ~ 8.67L" 



(107) 



and then we should search for a region of the parameters related to the singlet under the constraint. 

In our configuration, the modulus parameter of the Jacobi's elliptic function cn is determinec 
as a function of X = as in Eq. (73) and the relation between them is shown in Fig. [ll, 

where we understand that we have to take very small X for obtaining k = 1. We also refer to the 
complete elliptic integral of the first kind K[k], which is the function of the elliptic modulus k and 



whose value is equal to the quarter period of the Jacobi's elliptic function cn(y, k) in Fig. 11 This 
plot suggests that if we take the infinite period, whose case corresponds to l/cosh(?/, /c), we tune 
the value of k very close to one. 

Positions of divergent points of l/cn(x, fc), which correspond to zero points of cn(x, k), give us 
another important suggestion. The l/cn(a;, k) function gets divergent with a period of 2K[k\ and 
the range of [0, L] should not contain any such point. In the exact form of the VEV ($(y)) in 
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Figure 12: (a) the value of the scaled quarter period of the exact VEV in the function of X in 
the range of < X < 0.0001. (b) the value of the ratio, which is defined as the exact VEV form 



over its approximated form in Eq. (104). 



Eq. ( 75 ) , the position yd with divergence is evaluated as 



yd = yo + — 



1 



M(l +X) 



1/4 



K 



^(1 



mod — 



M(l + X) 



1/4 



K 



(108) 
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The second term of the right side in Eq. ( 108 ) means the quarter period of 



with X = X 

($(?/)) in the coordinate y. Considering the profile in Fig. |8| in the scaled coordinate y the position 
of ijd is preferred at one plus a few positive values. When we consider the property of the scaled 
quarter period in Fig. [12} as we have discussed before, we need to make the value of X approach 
to zero (but not exact zero) for 0{1) scaled period. In addition, we have to take account of the 



condition on yQ in Eq. (104), which is interpreted in the scaled coordinate y as 

Miy -yo)>l-^ 8.67{y -yo)>l in the range of < y < 1, 



(109) 



with the value of dimensionless M in Eq. (107). Here we observe that the case of positive yo is 
problematic (at least) around y = 0. Based on all the knowledge which we have obtained, we can 
find a set of parameters 



M = 8.67, 



yo 



-0.1, A = 0.001, |Q| = 0.001. 



(110) 



The validity of this choice is checked in Fig. [12] to calculate the ratio, which is defined as the exact 



VEV form over its approximated form in Eq. (104), and the difference is estimated as about fifteen 
percent at most. 

In the following, we check whether the EWSB occurs or not in our configuration. The pre- 



scription is written in the Section 4.2 and the two input parameters L± can be inversely calculated 



from the profile of the exact VEV through Eq. ( 76 ) as 

1 1 
^ = -6.07, ^ 



.69 



:iiii 
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and the value of Lmax in Eq- (78) is simultaneously fixed 1/L 



8.69 as the scaled values based 



on L. Agreeing with the previous naive discussion at the end of Section [4^ the condition 



LI 



(112) 



is fulfilled and then the EWSB is realized in our configuration for real. Here we comment on the 
two things. One is that the parameter A always appears in the cn function as the combination 
of \Q\\ and A in itself only affects on the overall normalization. This means that we can take 
more greater value in A with more smaller \Q\. The other is that the smallness of Uq and |Q| is 
not unnatural thing, because they are resultant values derived from the two input parameters 
whose dimensionless values are within (9(10). 

Owing to the above discussion, we are able to calculate Yukawa mass matrix elements in our 
model with the "elliptic VEV" ($(y)). We take care of the following two facts: 



with the parameters in Eqs. (110) and/or (111). 



We use the exact form of the VEV in Eq. (75) instead of the assumed warped form in Eq. (31 ) 



Because the Yukawa structure is introduced as the higher-dimensional operators in Eq. (98) 
the replacement is required 



The diagonalized Yukawa mass matrices take the forms of 

-M I diagonal = diag (2.47 MeV, 1.18GeV, 174GeV), 
-M^'^^ldiagonai = diag(3.94MeV, llOMeV, 4.19GeV), 

and the CKM matrix is given as 



CKM 



0.977 0.214 0.00448 
0.213 0.976 0.0475 
0.0145 0.0454 0.999 



(113) 



(114) 
(115) 



:il6) 



where we adopt the values in Eq. (Q for the 9 lengths {L^^\ L^, if, L^^\ L^^\ L^^\ if, if^.L^) 
and the 3 fermion bulk masses [Mq,Mii,Mj)). We can find only a small difference between the 
results and the previous results based on the assumed warped VEV. 

In contrast to the foregoing analysis, the Higgs dynamics is well described and the dimensionless 
coefficient part of ($(?/)), which is related to A in the assumed warped VEV, is calculable and the 
values of y^^'ly^^^ are as follows: 



= 0.0442, y^'^^ = 0.00369. 



(117) 



Here we mention that when we consider the situation where all the objects are localized, the 
effective length of each object is considered to be smaller than the total length L. If we take the 
effective length for overlap integrals showing down quarks as about L/3, which is the length of the 
integral for m^ii , the unnaturalness in y^'^^ is somewhat ameliorated as follows: 



3^^^^) ~ 0.00369 ■ (3)^ = 0.0332. 



118) 



The insertion of the relatively-large massive value v = 24 6 Ge V, compared to the quark masses 
except top quark's one, to the 5D Yukawa structure in Eq (113) is a cause of the smallness of y^"^^ 



and y^'''\ We do not go for more detailed discussion since we cannot avoid some ambiguities. 



26 



5 Summary and discussions 



We have presented a review on some properties of 5D fermion and vector fields in multi-interval 
system, where each interval is connected to the others through point interactions. By choosing 
BCs at the positions of point interactions suitably, the profiles of fermion zero modes get to be 
three-fold degenerated, localized, and mixed, which means all of the Yukawa structure in the SM 
can be realized. Combined with the ^/-position-dependent Higgs VEV with exponential warped 
form, the quark mass hierarchy and the structure of the CKM matrix are explained simultaneously 
almost only via geometry of extra dimension. One way to generate the warped VEV without gauge 
universality violation is to introduce both of the Higgs doublet with the ordinary Neumann BCs 
and the scalar singlet with the generalized BCs, which are coupled in higher- dimensional Yukawa 
terms. The ordinary Yukawa terms are prohibited by adding a discrete symmetry. The exact 
form of the scalar singlet VEV is written by the Jacobi's elliptic function and we have found that 
it becomes close to the exponential function in a region of parameters with almost (9(10) input 
parameters. To avoid violation in gauge universality, we should assume that the magnitude of the 



coefficient of the doublet-singlet mixing term in Eq. (99) is sufficiently small. 

Here we estimate the effect from KK mixing shortly. In our system, translational invariance 
along y is highly violated because of the existence of the point interactions and moreover KK-parity 
cannot be defined in behalf of lack of reflection symmetry. Consequently, the zero modes and KK 
modes of the fermions are mixed at the tree level and this affect the value of the mass eigenstates 
and the elements of the CKM matrix. Here we would like to consider the following form 



(0) 


















Mkk 






L 







Mkk 





-t(o)- 











+ h.c. 



:ii9) 



R 



where rrit, M^Kit^^KT^^^t^^^ are the top quark mass, KK scale, top zero mode, flrst top KK state 
of SU{2)]y doublet, that of SU{2)\y singlet, respectively. In the following, we only consider the 
mixing in the top quark sector because mixings in the other flve flavors are negligible because of 
smallness of their Yukawa couplings. In reality some deviation factors from the assumed ordinary 
UED-like form probably emerge in front of each matrix component originated from difference in 
overlap integrals and mixings between a zero mode and/or between zero modes and KK states of 
more than flrst level are also derived. But we ignore these issues for simplicity. 

The deviation ratio in the observed top mass, which is deflned as the mass eigenvalues over the 



reference value 172 GeV is calculated and the results is represented in Fig. 13 Then we conclude 
that we can ignore the level-mixing effect when the KK scale, which is deflned as vr/L in the 
ordinary UED context, is greater than 300 GeV up to accepting our naive assumptions. More 
detail analysis is out of the scope of this paper and we will discuss this issue in a future study. 
It is noted that in a "decoupling" case with a huge Mkk, we can always neglect the level mixing 
effect even if this possibility is not so interesting in collider physics point of view. 

The work done in this paper is considered as a flrst step for constructing a phenomenological 
model which explains the number of fermion generation, fermion mass hierarchies, and the structure 
of fermion mixing matrices simultaneously mainly based on geometry of an extra dimension. But 
lots of issues to be studied remain. 

The flrst issue is whether our mechanism works well in the lepton sector. As widely known, 
two mixing angles of neutrinos are large, and this suggests that the components of neutrino mass 
matrix are probably the same order of magnitude in contrast to the quark one. In addition, the 
expected neutrino masses are up to (sub-)eV order, where we find at least six orders of magnitude 
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Figure 13: The deviation ratio in the observed top mass, which is defined as the mass eigenvalues 
over the reference value 172 GeV. 

smaller than the value of the lightest charged lepton (electron). These differences are nontrivial 
and we will search for a configuration where the properties of the quarks and leptons are derived 
simultaneously via geometry. The second related issue is on the phase of CP violation in the 
CKM matrix. The existence of this phase with nonzero value has been established by B physics 
experiments, but within our present mechanism such a phase never occurs since all the zero mode 
functions are real. One considerable direction to overcome this problem is to introduce a complex 
phase through twisted boundary conditions. 

Some away from phenomenological issue, the Higgs with the generalized boundary conditions 
has a rich theoretical structure and there many topics wait for being unveiled. The phase structure 
of the scalar singlet is explored in Ref. [SO] but for non-singlet scalars some limited studies have 
been done. In non-Abelian gauge theory the boundary conditions can mix the gauge indices 
and possibilities are enlarged. Accordingly its phase structure gets to be much richer. Another 
interesting issue is the structure of quantum fiuctuation around the position-dependent elliptic 
VEV of the singlet scalar. Even in the case of the zero mode, the properties are highly nontrivial 
but this issue is mandatory when we discuss the signature of the scalar singlet at colliders. Moduli 
stabilization via Casimir energy is also an important topic for ensuring the stability of the system 
with the nontrivial VEV structure^] 

The physics in the system with point interactions and/or y-position-dependent scalar VEV just 
only starts to be dug out and we can find a lot of fascinating themes in both phenomenological 
and theoretical points of view. 
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Appendix 



A Estimating the orders of magnitude of doublet-singlet 
scalars mixing effect 

In this Appendix, we consider the minimization problem of the scalar potential of the doublet H 
and singlet $ scalars with the doublet-singlet mixing term in Eq. (99). In this analysis, we assume 
that the VEV of the singlet $ takes the effective form {4>{y)) in Eq. (31) and we concentrate on 
the part of 

j dyS^H\dyfH + M'^H^H-^{H^HY -C^^^H^Hy (120) 

After the replacement 



V2 



(121) 



we can identify the functional form £'[{h)] which we should minimize as follows: 



£'[{h)] 



dy< {dy {h)f - m" {h) 



(h)' + c {<pf {hy 



(122) 



Due to the Neumann BCs in Eq. (92), the form of the (position-dependent) VEV {h{y)) is fixed 



as 



-M' + l3o + J2(^^ cos ( 



n=l 



/irn 

It 



(123) 



with the coefficients /3o, (3n. Here the first term of the right-handed side of Eq. (123) corresponds 



to the solution with C = 0, whose value is equal to v/V2L with v = 246 GeV, and the remaining 
two terms show the deformation from it in the case of C 7^ 0. Under the assumption that the 
value of C is small, the potential is minimized with the forms of the coefficients 



/3o 



2M'L 2a 



(1 



-2aL\ 



((-1)" - 



(124) 
(125) 



within the second order of C. 

After we consider the suitable order estimation 



A' ~ L, M' ~ V, AVl ~ V, L-^ ~ Mkk, aL ~ C(l), 



(126) 
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where Mkk is a typical mass scale of the KK states, we can conclude that the orders of magnitude 
of the deviation from C = as 




C 

a 



a 



'KK 



(127) 



:i28) 



with the dimensionless values a = aL^^,C = CL. The results tell us that the value of /3„ is 
suppressed by the KK index n and Mkk, but on the contrary, that of f3o is not suppressed because 
of a = 0{1) in Eqs. (58) and (126) although /3o does not affect our conclusions since it merely 
shifts the constant expectation value of H. On the other hand, nonzero values of could cause a 
problem of gauge universality, so that C (Mkk) should be sufficiently small (large) in our model. 
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